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Starting from an initial state of thermal equilibrium, we derive an expression for the quantum 
fluctuation in the energy density during the inflationary epoch in terms of the mode functions 
j^ ' for the inflaton field. The effect of this particular initial state is not washed out in the final 

formula, contrary to what is usually believed. Numerically, however, the effect is completely 
negligible, validating the use of the two point function in the vacuum state. We also point out the 
1/^ ' requirement of conventional quantum field theory during infiation, that the quantum fiuctuation 

in a wavelength must be evaluated, at the latest, when the wavelength crosses the Hubble length, 
in contrast to the usual practice in the literature. 



PACS number(s): 98.80.Cq 



I. INTRODUCTION 



The most attractive aspect of inflationary models of the early universe m is their potential to predict the present 

day density inhomogeneity from first principles [g| . In these models it is possible to calculate quantum fluctuation 

^S^ ' in the energy density on the homogeneous background in a region within the causal horizon (given by the Hubble 

Ch, length) during the inflationary epoch. This fluctuation provides the initial spectrum of density perturbation. As 

Q,' the region inflates into the observed universe or bigger, its propagation through different eras can be followed till 

^ \ the present time by the equations of linear perturbation theory of classical gravity ||3|] . 

i-^ ' In this work we consider some points in the calculation of quantum fluctuation during inflation. The basic 

ingredient is the expectation value of the product of two scalar field operators at a time when considerable inflation 

has already taken place. It is generally believed that as the inflation proceeds, the effects of all scales associated 

with a particular initial state tend to be wiped out, retaining only the extremely high energies associated with 

quantum fluctuations in the vacuum. So the expectation value is evaluated for the homogeneous background, which 

corresponds to the vacuum state of quantum field theory. 

Here we investigate how far the density fiuctuation is actually independent of the initial condition prevailing at the 
beginning of infiation. For this purpose, we start with a definite initial state, namely that of thermal equilibrium. 
The first attempt in this direction was by Guth and Pi I?]. We discuss it here in a more general framework 0, H. 
The thermal propagator can be followed till the time when the quantum fluctuations are evaluated. 

The existence of an initial thermal equilibrium distribution of particles, at least for the high wave numbers needed 
for the calculation of density fluctuation, appears quite probable. Even if the collision rates among the particles 
are too small to produce such a state, there could be another mechanism at work. As pointed out by Weinberg |7| , 
the strong gravitational interaction at very early times would bring about thermal equilibrium, which, as we shall 
show, could be maintained at least till the beginning of inflation. 

The other point we discuss is the time at which quantum fluctuation must be evaluated. It relates to the 
applicability of quantum field theory in curved space-time. As emphasised by DeWitt B, conventional quantum 
field theory requires the mode functions to be oscillatory in time, allowing positive and negative frequencies to be 
identified. While on flat space-time such modes naturally arise for field theories describing physical particles, their 
existence is not guaranteed on space-times with non-zero curvature. The reason is that the curvature gives rise to a 
damping-like term in the equation of motion for the mode functions. In the inflationary period this makes a mode 
oscillatory or damped, according as the associated physical wavelength is smaller or bigger than the Hubble length. 
During this period physical wavelengths grow at a tremendous rate, while the Hubble length remains constant or 
approximately so. Thus even a wavelength lying initially deep inside the Hubble length would eventually go outside 
this length. So this time of exit marks the latest time at which we can evaluate the quantum fluctuation belonging 
to that particular wavelength. 

In the literature, however, quantum fluctuations are actually evaluated at a time, when the modes have evolved 
well outside the the Hubble length, so as to be frozen |9|. Of course, the complete problem of predicting the density 
fluctuation at the time of Hubble length reentry in a later radiation or matter dominated phase does involve the 



evolution of the fluctuation over a much longer period of time. But the question at hand is where the quantum 
fluctuation can be evaluated reliably. 

In Sect. II we review the derivation of the thermal scalar propagator in the early universe. In particular, we show 
the behaviour of mode functions as the wavelengths grow and discuss the validity of the assumed initial thermal 
equilibrium state. In Sect. Ill we write the formula for density fluctuation in terms of the mode functions and 
discuss its dependence on the initial condition. In Sect. IV we consider the simple, original model of extended 
inflation as an example ||l^]. Here we review the homogeneous classical solution for the inflaton fleld and set up 
quantum theory in this background. Finally, our concluding remarks are contained in Sect. V. 

II. FINITE TEMPERATURE SCALAR PROPAGATOR 

Consider a region of space in the early universe well within the causal horizon. It can then be taken to be 
homogeneous and isotropic, admitting the (spatially flat) line element, 

ds^ ^dt^ -a^{t)dx^ , (2.1) 

where the scale factor a{t) describes the expansion of the region of the universe. It constitutes the background 
space-time, which is perturbed by quantum fluctuations. The action for the scalar fleld in this space-time may be 
generally written as, 

S^ = ^fd^xdta^it) U^ - i(V0)' - fi^im^ ~ Xi{t)(t>' ~ X2m^ + ■■■]. (2.2) 

Here we have already shifted the scalar fleld by the classical, homogeneous field. The dots indicate interaction 
terms, if any, of the scalar field with other (gauge and matter) fields. 

We assume the different species of particles to be in thermal equilibrium around some initial time tg, which 
we conveniently take to be the time of transition of the radiation dominated phase to the inflationary phase. In 
particular, the scalar particles belonging to 4>{x) are also assumed to be in thermal equilibrium. (This assumption 
will be examined at the end of this section.) The density matrix is then given by 

where l//?o = T(to) is the temperature at time tg. The explicit time dependence of the Hamiltonian Ti.{t) arises 
from that of the scale factor and the homogeneous classical field. Note that the density matrix is constant in 
the Heisenberg representation. Thus once the system is in the thermal equilibrium state, the thermal propagator 
continues to hold even when the system deviates from this state. 

To describe the time evolution of the system, it is most convenient to use the real time formulation of thermal 
field theory [Q. In the context of the early universe, the time path C in the action integral consists of three 
segments as shown in Fig.l [||. The points on it may be labelled by a complex parameter r such that 

t, on Ci and C2 ^^ 4) 

to — it, on C3 . ^ ' ' 

The action in the path integral corresponding to the segments Ci and C2 is in Minkowski space, while it is 
Euclidean on C3. It should be noted, however, that the scale factor is not continued to Euclidean space on C3 : 
since the Hamiltonian for the segment is 7i(to), the scale factor remains fixed at tg. 

We now review the derivation of the thermal propagator |q| , [p| . After a partial integration, the quadratic terms 
in 5*0 becomes [[12| , 

5*0 = — / dr d^x(j)D(f) + boundary terms, (2.5) 



where 



— -3ti^+c.2 , ^2(^)^_ 1 v2+^2(^)^ reCuC2 



We use the abbreviations, ao — a(io), i-^Q ~ "-^(^o), Tq — T{to). For the boundary terms to vanish, it is necessary 
that both (j) and -^ match at the joining of the segments Ci and C2, of C2 and C3 as well as at the free ends of C3 
and Ci. 

The (time ordered) thermal propagator TrpT(j){x)(j){x ) will be denoted by Gi3{x,t]x ,t ) or < T(j){x)(j){x ) >. 
It satisfies 



DGfsix, T-x ,t)^ -i5^{x - X )6{t - t ) , (2.7) 

with boundary conditions following from the matching of and -^ mentioned above. For the spatial Fourier 
transform of the propagator, defined by 

^^e*-G^(fc,r,r), (2.8) 

it reduces to 

DG0{k;T,T')=-i6{T~r'), (2.9) 

where — V^ appearing in the expressions for ^^^(t) and uJq is to be replaced now by k^. 

To construct the propagator we first find the mode functions. On the contour C3 they satisfy 

^+,oiy^iT) = 0, (2.10) 

givmg 

/i±(t) = -=L=eT'"°^ T = t„-tt eCs. (2.11) 

V2woao 

The normalization satisfies the Wronskian condition, }i^{T)h^{T) — h^ {T)h^ (t) = —i/a^. The mode functions on 
the real segments Ci and C2 are the solutions of 

iL+3^A+^2(^)yg±(^)^0, r = ieCi,2, (2.12) 

ar'^ adr J 

with normalisation fixed again by the Wronskian condition, (j^ {T)g^ {r) — g^ {T)g^ (t) = —1/0^(7). To see the 
nature of these solutions, we put 

5±(r)=a-3/2-±(^)^ (2.13) 



where g satisfies 



with 



|^+cD2(r))g±(r)=0, (2.14) 



9.x k"^ 9 9/0 2 • \ , , ait) 



For a power law behaviour of the scale factor, a{t) ~ t^, it becomes 



-.2/j.\ ''' , ..2 9 (1 2 \ rr2 



^^t)^-,^,^--[l--jH^. (2.16) 

It is now simple to identify the modes, which belong to conventional quantum field theory. The magnitude of 
yLt(t) is usually small compared to H{t). Thus in the radiation dominated phase [p = 2)1 ^'^ is positive for all values 
of k and it may be possible to define oscillatory modes belonging to positive and negative frequencies, at least in a 
quasi-static way. We show below that this is indeed the case around the time ig, when we can solve (2.14) in the 
JWKB approximation to get |13|, 

I . , 1 =Fi / dr lu(t ) 

ff±(r) = -===e^ J'o >, r^to. (2.17) 

We thus have a valid quantum field theory around the time io- 

But in the inflationary phase (p 3> 1), a mode is oscillatory only if its physical wavelength 2TTa{t)/k is small 
compared to the Hubble length H^^{t). As inflation progresses, the scale factor increases enormously, while H{t) 
is approximately constant. Thus more and more modes go out of Hubble length and behave as damped waves, 
having no interpretation in quantum field theory. 

The solutions g^ir) and h^{T) may now be joined to form the functions /^(t) on the entire contour C, 



By definition, /^(t) obey the continuity conditions relating the segments Ci and C2. Using Eqs (2.11) and (2.17), 
we see that the conditions connecting C2 and C3 are also well satisfied if H{to) is small compared to k/a{to) p| . 
A particular solution to (2.9) may now be written as 

Goik ■.t,t') = f+{T)f-{T')e,{T - r') + f+{T')r{r)ec{T' - r), (2.19) 

where 9c is a step function on the contour. It satisfies the continuity conditions at the junctions of segments Ci 
and C2 as well as of C2 and C3, because /^(t) does it. To satisfy the remaining (thermal) continuity condition at 
the ends of Ci and C3, we add to it the most general solution of the homogeneous equation, 

2 
Gp{k;T,T') = Go{k,T,T') + ^ r{T)A^^f{T') . (2.20) 

The superscript (±) on the mode functions are replaced temporarily by 1 and 2 to use matrix notation. The 2x2 
constant coefficient matrix A is uniquely determined by the thermal conditions |1^]. We get 

Gp{k,T, t) = /+(r)/-(r'){0,(T - r') + n{u:^)} + /-(r)/+(T'){0,(r' - r) + n(c.o)} , (2.21) 

where n(wo) is the bosonic distribution function 

n(wo) = (e'^°'^° - 1)-^ (2.22) 

For tree level calculations we need the Green's function only on the real axis Ci. Writing henceforth t = t, this 
is given by 

{cj,[x,t)cj,{x,t))^j-^e^^-^^-'^\l + n{uJo))9+{t)9-{t), t>t, (2.23) 

where the mode functions g^{t) are solutions of (2.12). 

We now come back to the assumption of the initial thermal equilibrium state. Such an initial state can be 
ensured in an expanding universe if collisions among particles occur at a rate faster than the expansion rate of the 
universe. While this condition holds for species interacting through (relatively large) gauge coupling, it may not 
hold for particles of the inflaton field, which is a gauge singlet and has weak self-interaction. We discuss below the 
other mechanism, mentioned in the introduction, which could give rise to thermal equilibrium around the time to. 

At the Planck time tp, the strong gravitational interaction brings all species into thermal equilibrium M. If the 
system is quantised at this time in a cubic volume with sides small compared to the Hubble length, the longest 
wavelength will be well inside this length, i.e. 

> TTH{tp), (2.24) 



a{tp 
even for the smallest wavenumber. Then eq (2.16) simplifies to 



^{tp) = -^y (2.25) 



to a good approximation and the density distribution becomes 



"('^(*^)) = efe/a(t.)T(.p) - 1 ■ (2.26) 



We now bring the inequality (2.24) to the time to 



k ^ ^E^^Hito) ^^nHito), (2.27) 



a{to) a{to) H{to) ' ' To 

where mp is the Planck mass. In the last step we have used the radiation dominated solution for a{t). (See eqns. 
(4.6-7) below.) The temperature Tq is given by the grand unification scale, Tq ~ 10^^ GeV , so that mp/To ~ 10^. 
Thus the relation (2.25) at time tp continues to hold throughout the radiation dominated phase; in fact, it becomes 
more and more accurate as t increases from tp. Clearly the equilibrium distribution (2.26) established at time tp 
is well maintained at least till the time tg. 



III. DENSITY FLUCTUATION FORMULA 



The density inhomogeneity (at time t) is measured by the mean square fluctuation in the density function p{x, t) 

spY _ / f p{x,t) - p{t)Y\ 



pj \ V pit) 

where (• • ■)x denotes averaging over space and p is the homogeneous background density, p = {p{x,t))x- In the 
inflationary scenario, this fluctuation in the early universe is supposed to arise from quantum fluctuation. We may 
calculate the latter by evaluating an expression similar to (3.1), with p{x, t) replaced by the corresponding operator 
p{x, t) and the averaging process by taking the expectation value in an appropriate state. There is, however, a 
technical problem with this quantum version, as it involves the product of p{x, t) with itself at the same space-time 
point, which is not defined in quantum field theory. The problem may be avoided by taking the smeared density 
function [Q, 

pi{x,t) = N j d'ye-y"/'-'" p(x + v,t), (3.2) 

where I is an arbitrary smearing length and N an irrelevant normalization factor. The fiuctuation in pi is given by 

'5pi\'^ _ / f pi{x,t) - pi{t)\^\ 



Pi J c \ V Piit) 

where the subscript c stands for classical. The corresponding quantum fiuctuation, denoted by the subscript q, is 
now well defined, 



hiS^ ( Pi{.x,t) - p{ty 



Pi J a \\ P(t) 



(3.4) 



where (• • •) stands for the expectation value in the initial thermal state defined by eqn (2.3). 
To treat perturbation on different length scales, one writes 

p{x,t)^p{t){l + S{x,t)), (3.5) 

and Fourier analyzes the so-called density contrast, S{x, t), 

5(f,t) = ^^4(0e''^ (3.6) 

where 1^ is a volume within the Hubble length to begin with. In the limit of large volume, (3.3) becomes 



5pi \ f (Pk 



Mt^e-' ' . (3.7) 



PI J, J (2^)3' 
The energy density operator p(x, t) is the time-time component of energy momentum tensor, 

T^. = a^^a.* - 9^. Q.g'"'a«*9/3* - 1^(*)') , (3.8) 

for the full scalar field 'i'{x). Here the potential function V{'i) depends on the model considered. We shift \l/(x) 
by the homogeneous classical field ip{x), 

*(x,t) = V(t) + 0(^,i), (3.9) 

such that for the quantum field < (j){x) > — 0. The p{x) is given by 

pix)=pit) + Uix), (3.10) 

where 

m^li^^ + vw, (3.11) 



and 

U{x)^r{t)<P{x)+s{t)^{x), (3.12) 

to first order in (/'(a;). The coefRcients r{t), s{t) in (3.12) depend on the classical field and other parameters in the 
potential V^(^). Terms in U, which are of higher order in would give loop contribution to the density fluctuation 
and are neglected. Using (3.10), (3.12) and (2.23), the expectation value in (3.4) may be evaluated to give. 



spiV 1 



PI J a P'it) 



d^xd^ye-'-'^'+y'^^^^ < U{x,t)U{y,t) > 



pHt)J (2^)^ 

Comparing eqs (3.7) and (3.13) we get the desired result, 

1 



^ ^ '^''' (1 + n{ujo)) I rit)gkit) + s(i)5fc(t) P e-'^''' . (3.13) 



Skit) 1"= ^(1 + nitoo)) I r{t)gk{t) + s{t)gk{t) {' . (3.14) 

P 



Following our discussion in Sect. II, we evaluate it at time th, when the wavelength crosses the Hubble length, 

k 



a{th) 



2t: H{th 



The treatment of the evolution of density perturbation outside the Hubble length is standard M& , |I^ . At the 
time of its reentry within the Hubble radius of the post-inflationary radiation dominated epoch, it is given by [p_8|. 



Sp\ _ y^k^{l + n{LJn)) \ rgk{th) + sg ujih) \ ,^ ^^. 



pJh 3V27:{p + p)t^ 

where p is the homogeneous pressure, p — ^tp'^ — V{tp). 

It is simple to estimate n(ct'o) in the range of fc/oo, which is of interest. We write 

ao a{tp) a{te) a{to) 

From the time te when the inflation ends till the present time tp, the universe expands adiabatically, so that 
a{tp)/a{te) — To/Tp. The other ratio a(te)/a(io) = Z gives the magnitude of inflation. We thus get 

k 2Tr Z I Z , , 

(3.17) 



aoTo X{tp) Tp Xmpc 10^^ 

where Tp = 2.7K and Ajv/pc is A(<p) expressed in AIpc. The wavelengths of interest stretch over the range 1 < 
^Mpc < iC. To solve the problems of the standard cosmology we need Z > 10^^. But actually Z exceeds this limit 
by many orders of magnitude in most of the models of inflation. Thus k/aoTQ is large in these models and we may 
set n{ujo) = in the expression (3.14). 

We thus see that although the initial thermal equilibrium state does produce a factor in the expression for the 
density inhomogeneity, its magnitude turns out to be unity, justifying the use of zero temperature propagator for 
its evaluation. Nevertheless it is important to know the initial state, as there are other quantities, such as the 
duration of inflation, which may depend sensitively on it. 

IV. EXAMPLE OF EXTENDED INFLATION 

In this section we consider the well-studied, original model of extended inflation |1Q|. We flrst summarise, 
for the sake of completeness, the classical inflationary solution for the scalar field, which joins smoothly to the 
radiation dominated solutions before and after the inflation. We then write the Lagrangian for the quantum field 
in the background of this homogeneous field. The mode functions and the (vacuum) propagator has already been 
evaluated in a recent work |jl8] . Here we only discuss the possibility of establishing thermal equilibrium through 
collisions. 

The model is based on the Brans- Dicke (BD) theory of gravity given by the action M^ , 



The time dependence of the BD field $ makes the effective gravitational 'constant' vary with time, (jj is a dimen- 
sionless parameter. Cmatter represents the contribution of all other fields including the inflaton field ct, 

Cmatter = ^g^^'d^ad.a - V{a) + • • • . (4.2) 

This simple model is not realistic, however. The bubble nucleation problem can be solved for lo < 25, while 
astrophysical observation constrains it to lj > 500. Our purpose here is to examine in this toy model the possibility 
of realising thermal equilibrium as the initial state for the inflationary epoch. In the following we take a; > 25, say. 
The equations satisfied by the classical fields, viz. the scale factor a{t) and the homogeneous part (p{t) of the BD 
field $(x) are |l§, 

(4.3) 



^ + 5H^=^^{p-3p) , H^ 


a 
a 


2 Svrp u fifiV fip\ 
3ip 6 \ipj \ipj 





(4.4) 

where the energy density p and the pressure p are generated by Cmatter ■ If the masses are small compared to the 
temperature, they are 

P=^NT^ + M\ p=^NT^~M\ (4.5) 

where Af^ — V{0) is the false vacuum energy density and N counts the total number of effective degrees of freedom. 
In the pre-inflationary radiation dominated epoch, (p{t) is constant [po[ , 

^(t) = ifiito) , ait) ex Vt, (4.6) 

giving 

(4.7) 



(4.8) 











Th^ 


/45</p(to) 
V 167r3iV ' 
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epoch the solutions 


may be 


written as 













Vit) 


^ipito){l + B{t- 


to)r, 




with 






B^ito) = 


M2 


,-^^ 








3w + 5)(2w 


+ 3) 



and 



a{t)^aito){l + B{t-to)r+'/' . (4.9) 



The two sets of solutions join at io: when we have approximately 

NT^ ~ M^, (4.10) 



7r2 



30 
giving To ~ M and 



.2, _ . / 45</?(to 
167r3iV 



To%-\h£;^- (4-11) 



Equating the values of -ff (to) for the two solutions, we get another relation 

^^(.- + i)B, («2) 

It is no new constraint, coinciding with (4.11) to leading order in uj. 

We now assume as usual that the inflationary period ends instantly at time t^ in radiation domination again at 
a temperature T ^ Tq. Then ip{t) ceases to vary appreciably, so that we have 



(p{t)^^{Q=mj,, t>te, (4.13) 

where mp = G^^^^ is the present value of the Planck mass. Using (4.8) and (4.9) we can relate ip{to) to the amount 
of inflation Z, 

^(to) = n4z-4/(2-+i), ^=^. (4.14) 

a(to) 

To set up a semi-classical quantum theory in the background of the above solutions, one faces the problem of 
dealing with the non-standard form of the action (4.1). This problem may be avoided by going over to a new ( 
Einstein ) frame from the present ( Jordan ) frame by a conformal transformation of the metric. Denoting quantities 
in the new frame by a bar, the required transformation is [pi| , 

g^,{x,t)^n-\t)g^,{x,t), (4.15) 

where 

r)2(i)^^. (4.16) 

^ ^ $(t) ^ ^ 

One also introduces a field $ to replace $, 






* = xln(— ), x=\\^^fnp. (4.17) 



which brings the kinetic term for $ in the canonical form. Assuming u — constant in the inflationary phase, the 
new action becomes, 

S = Jd'x^i^j^^ + ^g^^d^^d,^ + ^(*)| , (4.18) 

where V{'^) = M^e~^^^^. In the Einstein frame $ plays the role of the inflaton field. 

The new metric may be brought back to the Robertson- Walker form by a redefinition of the time coordinate 

dp - d^t)dx^ = n~^{t){dt^ ^ a\t)dx^} , (4.19) 

giving 

di=n-\t)dt, a{t)^Q-\t)a{t). (4.20) 

Using (4.16) and (4.20), the inflationary solutions (4.8-9) in the Jordan frame may be recast in the Einstein frame. 
Integrating the first equation in (4.20) and choosing the constant of integration properly, we get |22[| , 

l + C{i-io) = {l + B{t-to)r, C=^^^. (4.21) 

V¥'(io) 

Then the new scale factor is obtained from the second equation in (4.20), 



a{t) = a{io){l + Cii- io)}(2-+3)/4, a(io) = a(to)^^^ • (4.22) 

rap 

The classical homogeneous part of \l/ is obtained from (4.8) and (4.17), 

m = Hto) + Xln{l + Citato)}, V(io) = Xln (^) ■ (4.23) 

These solutions could, of course, have been found by solving the equations obtained by varying the action (4.18) 
in the Einstein frame. 

The free part of the action (4.18) is now in the standard form for developing quantum field theory. We introduce 
the quantum field (j){x) by decomposing ^ as, 

*(f,i)=V(i)+ 0(^,0- (4.24) 

The exponential potential V{'^) as such has no place among renormalizable quantum field theories. But we can 
treat it as an effective potential and retain terms up to fourth power in cj). The action for cj) then becomes 




where 



S = / d'xdia^ {^[^] - ^(V0)' - ^m'0' - Ai0=^ - A204 I , (4.25) 



l^H^--H^it), (4.26) 



Ai(i-) 
A2(i) 



2r 



167ri72(f) 



The magnitudes of /i^, Ai and A2 depend on l/(p{to), the effective gravitational constant in the pre-inflationary, 
radiation dominated phase. The relation (4.14) is not useful in determining <p(to) as Z itself depends on the details 
of the models of inflation. However a lower bound on (p{to) may be obtained from an earlier requirement that 
fi{to) < To, giving 



/ N /327r ,^ 
<p(io) > \ Mmp. 

V IjJ 

2 / \^ 

From the corresponding bounds on Ai and A2, it is simple to find that the collision rate, Tcoii < (|^) ( ;^ ) M , 

which is hopelessly small compared to the expansion rate, T^^p = H{to). Even though collisions are totally 
ineffective to maintain thermal equilibrium, wc have shown in Sect. II that gravitational interaction at very early 
times would ensure equilibrium distribution of the (p particles around the time ig- 

There are several earlier calculations of density fluctuation in this model [E2] , p3] . All are based on the result for 
the vacuum propagator for the scalar field, which we justify in Sect. III. However earlier works generally evaluate 
the quantum fluctuation outside the Hubble length, which following our discussion in Sect. II, cannot be interpreted 
properly in quantum field theory. In a recent work |lq | , we evaluate this density fluctuation at the time of Hubble 
length exit using the formula (3.15). We find a result, which is an order of magnitude bigger than the earlier ones. 

V. CONCLUSION 

In the present work we assume the inflationary epoch to begin in a state of thermal equilibrium and study its 
effect on the quantum fluctuation in the energy density calculated during this epoch. This initial state including the 
scalar particles appears quite likely even if their self -interaction is too feeble to ensure it. We show that the thermal 
equilibrium established at very early times through the-then strong gravitational interaction would be maintained 
till the beginning of inflation. By evaluating the scalar field propagator with thermal boundary conditions, we find 
a result for the density fluctuation, which differs from the one calculated with the vacuum propagator by the factor 
{1 + (e'"'"/'^"^" — 1)^^}. Clearly the factor does not go to unity as time passes but is a constant depending on the 
physical wave number and the temperature referred to the initial time io- 

It turns out, however, that for wave numbers of interest in the present universe, this factor is unity in models 
where the amount of inflation exceeds by many orders of magnitude the minimal amount required to solve the 
problems of standard cosmology. Thus numerically the calculation of fluctuation in the vacuum state is justified. 

We also point out that the conventional quantum field theory applies on curved space-time as long as the modes 
oscillate. This requires that we evaluate the quantum fluctuation, at the latest, when the corresponding wavelength 
crosses the Hubble length. Previous works 0, however, evaluate it as a rule for wavelengths well outside this length, 
where the modes freeze. As we have shown in a recent work |1^ , this difference in the calculation leads to an increase 
in the result by a factor of five for the model of extended inflation. 

Finally we note that as long as the modes are within the Hubble length, they retain a thermal equilibrium 
distribution. Thus although the initial state of thermal equilibrium comprising all modes is not maintained during 
inflation, the modes relevant for the calculation of density fluctuation are those still in an equilibrium distribution. 
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Fig. 1. Time path of real time thermal field theory 



